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1.  INTRODUCTION 

The  primary  objective  of  the  research  conducted  under  this  grant  is  to  advance  funda¬ 
mental  understanding  and  prediction  capabilities  of  reacting  flows  relevant  to  Air  Force 
needs  in  high-speed  airbreathing  propulsion.  The  dynamical  coupling  among  the  tur¬ 
bulence,  thermodynamic  variables  and  fuel  chemistry  is  not  well-understood.  Advanced 
simulations  of  prototypical  shear  flows  and  associated  analytical  tools  are  used  here  to 
understand  and  to  model  important  aspects  of  this  coupling. 

The  computational  tools  are  primarily  direct  numerical  simulation  (DNS)  and  large 
eddy  simulation  (LES).  In  DNS,  the  conservation  equations  for  a  compressible,  multi¬ 
species,  reacting  mixture  are  numerically  solved  to  capture  all  dynamically  important 
scales  of  turbulent  motion  without  resort  to  any  turbulence  closure.  In  LES,  filtered 
equations  that  represent  the  behavior  of  the  ‘large  eddies’  are  simulated  and  the  influ¬ 
ence  of  the  unresolved  ‘small  eddies’  modeled  with  an  appropriate  subgrid-scale  (SGS) 
model.  Both  DNS  and  LES  involve  the  numerical  solution  of  unsteady,  three-dimensional 
forms  of  the  governing  equations.  LES  requires  fewer  computational  grid  points  than 
DNS  and,  provided  that  accurate  SGS  models  are  developed,  can  be  used  for  high- 
Reynolds  number  flows.  DNS  has  the  advantage  of  giving  results  that  do  not  depend 
on  possible  flaws  in  turbulence  closures  and,  therefore,  is  best  suited  for  developing  new 
insights  and  uncovering  new  phenomena. 

Th  research  accomplished  as  part  of  this  grant  has  resulted  in  19  publications  and 
two  Ph.D.  theses. 

2.  THE  COMPRESSIBLE,  REACTING  SHEAR  LAYER 

DNS  has  been  used  to  study  in  detail  the  compressible,  reacting  shear  layer  and 
to  discriminate  among  compressibility  effects  due  to  high-speed,  variable  density  and 
exothermicity.  The  flow  consists  of  hydrocarbon  fuel  and  air  streams  with  a  relative 
shear  which  mix  and  react.  Unlike  previous  simulations,  large  computational  domains 
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Figure  1:  Dependence  of  the  shear  layer 
growth  rate  on  the  convective  Mach 
number,  Mc.  The  growth  of  the  shear 
layer  thickness,  8,  is  normalized  by  its 
incompressible  value. 


Figure  2:  The  dependence  of  turbulent 
kinetic  energy  production,  P,  and  dissi¬ 
pation,  e  on  Mc.  The  production  and 
dissipation  are  normalized  using  AU 
and  6g. 


are  used  along  with  broadband  initial  fluctuations  to  allow  evolution  to  the  self-similar 
state.  Indeed,  it  is  found  that  qualitative  differences  in  the  turbulence  structure  and 
dynamics  exist  between  the  self-similar  shear  layer  and  the  earlier-time,  developing  shear 
layer.  Grids  with  up  to  13  million  points  are  used.  The  final  Reynolds  numbers,  based  in 
vorticitv  thickness,  are  large  enough  for  turbulent  flow,  reaching  values  as  large  as  12, 000. 
High-order  discretizations  are  used  to  ensure  numerical  accuracy  with  fourth-  and  sixth- 
order  compact  schemes  for  the  spatial  derivatives  and  a  third-order  Runge  Kutta  scheme 
for  time  advancement.  The  integral  length  scales  are  sufficiently  small  compared  to  the 
dimensions  of  the  computational  box  to  have  accurate  large-scale  representation,  while 
the  grid  step  is  sufficiently  small  to  resolve  the  small  scales. 

Highlights  of  our  results  are  summarized  below  with  further  details  available  in  arti¬ 
cles  [1,  2,  3,  4,  5,  6,  7,  8]  and  a  Ph.D.  thesis  [9]  resulting  from  this  work. 

2.1  The  high-speed  effect 

The  convective  Mach  number,  Mc  =  (U\  —  U2)/ (ci  +  C2),  is  often  used  to  characterize 
high-speed  shear  layers.  Here  U\,  cx  denote  the  velocity  and  sound  speed  in  the  upper 
free  stream,  while  U2,  c2  refer  to  the  lower  stream.  The  dramatic  reduction  of  shear 
layer  growth  rate  with  increasing  Mc  is  well-known.  As  shown  in  Figure  1,  our  DNS 
results  obtained  with  equal-density  streams  are  in  close  agreement  with  the  Langley 
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experimental  curve  and  other  data  sets  which  use  air-air  streams.  The  experimental 
data  points  with  reduced  growth  rates  with  respect  to  the  Langley  curve  correspond  to 
heterogeneous  gases  with  unequal  density. 

With  DNS,  we  are  able  to  uncover  a  mechanism  which  accounts  for  the  reduction 
in  the  shear  layer  growth  rate  with  increasing  Mach  number.  The  decreased  level  of 
turbulent  kinetic  energy  production  (see  Fig.  2)  is  found  to  be  responsible  for  the  reduc¬ 
tion  in  question  in  agreement  with  the  conclusion  drawn  by  [10]  from  DNS  of  uniformly 
sheared  flow.  The  normalized  production  has  been  also  found  to  decrease  with  Mc  in 
an  independent  study  of  the  shear  layer[ll]  and  in  the  annular  shear  layer  [12].  Fig.  2 
shows  that,  while  the  normalized  turbulent  production,  P,  decreases  with  increasing 
Mach  number,  the  turbulent  dissipation,  e,  does  not  decrease  proportionally,  resulting 
in  a  net  decrease  of  turbulent  kinetic  energy. 

Further  examination  of  the  Reynolds  stress  transport  equations  using  the  DNS 
database  shows  that  a  reduction  in  the  pressure-strain  term  reduces  the  energy  transfer 
from  the  streamwise  component  to  the  other  components,  thus  decreasing  the  Reynolds 
shear  stress  and,  thereby,  the  turbulent  production,  P .  The  pressure-strain  term  is  re¬ 
sponsible  for  transfer  of  energy  among  the  three  velocity  components  and  is  defined  by 
n,j  =  rfidu'Jdxj  +  du'j/dxi),  where  p‘  is  the  pressure  fluctuation  and  the  velocity 
fluctuation.  Fig.  3  shows  that  all  components  of  the  pressure-strain  tensor,  II decrease 
with  increasing  Mc.  The  normalized  pressure  fluctuations  (not  shown)  also  decrease 
with  increasing  Mc. 

An  analysis  to  uncover  the  mechanism  underlying  the  reduction  in  the  pressure-strain 
term  has  been  performed.  In  strictly  incompressible  flow,  the  instantaneous  pressure 
has  ‘infinite’  signal  speed  and  satisfies  a  Poisson  equation  whose  source  is  related  to  the 
velocity  gradients.  In  the  case  of  compressible  flow,  pressure  fluctuations  travel  with 
the  speed  of  sound,  cq.  The  consequence  of  a  ‘finite’  signal  speed  on  the  pressure-strain 
term  is  studied  analytically. 

The  following  evolution  equation  for  the  pressure  fluctuation  p  , 


1  D2p  d2p  _  d  /  \ ' 

c2  Dt 2  dxidx  i  dxidxj  v^1  j) 


(1) 


can  be  derived  by  taking  the  divergence  of  the  momentum  equation,  and  assuming 
that  the  isentropic  relationship,  Dp/Dt  =  c20Dp/Dt ,  applies  and  that  viscous  terms  are 
negligible.  Such  assumptions  are  reasonable  for  high-Reynolds  number  flows  away  from 
shocks  and  solid  boundaries.  Eq.  (1)  is  now  specialized  to  the  center  of  the  shear  layer 
where  the  mean  velocity  is  zero,  du\/dx2  ~  constant,  and  an  analysis,  assuming  locally- 
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Figure  3:  Dependence  of  components  of 
the  pressure-strain  tensor,  lip,  normal¬ 
ized  by  corresponding  incompressible  val¬ 
ues,  on  Mc.  Subscripts  1  and  2  denote 
streamwise  and  cross-stream  directions, 
respectively 


Figure  4:  The  performance  of  the  new 
pressure-strain  model  in  predicting  the 
compressibility  effect  on  the  shear  layer 
growth  rate 


homogeneous  turbulence,  is  performed  that  involves  two-point,  two-time  correlations. 
Details  of  the  analysis  are  given  in  [1].  The  result,  for  low  Mach  number,  is 

=  1  -  K0Af!  (2) 

ilij 

which  shows  that  rip  exhibits  monotone  decrease  with  Mc  in  agreement  with  the  DNS 
results  shown  in  Fig.  3. 

The  wave  equation,  Eq.  (1),  for  the  pressure  is  a  useful  approximation  is  consistent 
with  the  energy  equation,  if  heat  conduction  and  viscous  dissipation  terms  are  assumed 
negligible  and  the  isentropic  relation  is  assumed.  However,  since  these  assumptions 
are  not  physically  supported  when  Mc  is  not  small,  at  high  Mach  number,  the  general 
thermodynamic  relation  must  be  used, 

Dp  _  (dp\  DR  +  (^R\  —  (3) 

Dt.  \dp)sDt  \dS)PDt 

where  the  first  term  on  the  r.h.s.  is  the  ‘acoustic’  contribution  and  the  second  term  is 
the  ‘entropy’  contribution.  Indeed,  by  taking  the  time  derivative  of  Eq.  (3)  and  using 
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the  divergence  of  the  momentum  equation,  a  modified  wave  equation  for  the  pressure 
fluctuation  results,  which  involves  entropy  fluctuations  and  is  valid  for  all  Mach  numbers. 
Analysis  of  the  modified  wave  equation  shows  that  Ily  — >  const ,  for  large  Mc.  Thus, 
analysis  provides  functional  forms  for  small  and  large  Mach  numbers. 

A  predictive  model  for  the  pressure-strain  correlation  which  has  the  analytically 
predicted  behavior  at  low  and  high  Mach  number  is  obtained  by  using  the  ‘experimental 
values  of  nu  deduced  from  the  Langley  growth  rate  curve  to  determine  the  coefficients 
in  the  asvmptotic  solutions  and  using  a  exponential  fitting  function  to  approximate  the 
difference  at  intermediate  Mach  numbers.  The  new  pressure-strain  model, 


ni7  „  N1  +  aMc2e-(2Mc'1/2): 

nT  =  (1"c) - iTkA/j - 


+  C 


(4) 


for  the  compressibility  effect  on  the  pressure-strain  term  can  provide  the  dependence  of 
growth  rate  on  Mc  as  shown  in  Fig.  4. 

Based  on  analysis  and  DNS  observations,  the  following  mechanism  is  proposed  to 
explain  the  fundamental  property  of  compressible  turbulence  that  results  in  the  reduced 
growth  rate  of  the  shear  layer.  The  finite  speed  of  sound  implies  that  the  temporal 
decorrelation  (over  and  above  spatial  decorrelation)  between  adjacent  points  becomes 
additionally  important  with  increasing  Mc.  Consequently,  the  pressure-strain  term  de¬ 
creases  and  the  turbulence  redistribution  necessary  to  maintain  shear-driven  turbulence 
decreases  so  as  to  reduce  the  thickness  growth  rate. 

The  anisotropy  of  the  Reynolds  stresses,  =  pu-u'  /p,  is  an  important  characteristic 
of  the  velocity  fluctuations  used  in  advanced  modeling  of  turbulent  flows.  The  tensor, 
bt}  =  Rtj/2K  -  5ij/Z ,  where  K  is  the  turbulent  kinetic  energy  is  used  to  characterize 
Reynolds  stress  anisotropy.  The  effect  of  Mc  on  the  anisotropy  is  controversial,  with  some 
investigations  reporting  a  substantial  increase  of  the  anisotropy  in  the  high-speed  regime 
and  others  reporting  negligible  changes.  Two  important  conclusions  are  drawn  from 
our  DNS.  First,  during  its  initial  evolution,  the  anisotropy  tensor  is  strongly  affected 
by  Mc  such  that  larger  magnitudes  of  the  diagonal  components  of  bij  are  measured 
for  increasing  A/c.  Second,  after  sufficiently  long  time,  the  values  of  the  anisotropy 
tensor  approach  asymptotic  values  only  weakly  dependent  on  Mc.  The  peak  values 
of  turbulence  intensities  and  shear  stress  at  the  center  of  the  shear  layer  can  also  be 
considered  and  indeed  are  preferable  in  comparing  DNS  with  experimental  data.  The 
new  results  obtained  by  [13]  who  measure  all  three  turbulence  intensities  as  well  as  the 
Reynolds  shear  stress,  and  the  older  data  of  [14]  are  in  good  agreement  with  the  present 
DNS. 


5 


2.2  Effect  of  density  ratio 

The  convective  Mach  number  has  often  been  assumed  to  be  the  sole  parameter  that 
determines  compressibility  effects  in  the  shear  layer.  However,  there  is  significant  scatter 
in  Fig.  1  among  the  air-air  data  (Langley  experimental  curve,  Elliot  and  Samimy  (1990), 
Debisschop  and  Bonnet  (1993),  Chambres,  Barre  and  Bonnet  (1997))  and  the  other  data 
sets  which  used  gases  with  different  densities.  DNS  is  performed  to  investigate  potential 
effects  of  the  density  ratio.  The  ratio  of  the  freestream  densities,  s  =  fc/pi,  is  varied- 
while  the  average  density  po  =  (pi  T  P2 )  /  2  =  0.7  constant.  The  convective  Mach  number 
is  also  held  fixed  at  Mc  =  0.7.  Values  of  s  =  1,  4  and  8  are  examined. 

The  momentum  thickness  defined  by 

1  _  WA U  _  v  ,  /r\ 

60  =  — —  2  I  p(~z - Ui)(— +  Ui)dx2  .  (5) 

p0Au2  </_0O  2  2 

and  the  vorticity  thickness  defined  by 

5W  =  Au/(dui/dx2)  max  (6) 


are  popular  measures  of  the  shear  layer  thickness.  Here  AU  is  the  velocity  difference 
across  the  shear  layer  and  x2  is  the  cross-stream  coordinate.  It  is  important  to  distinguish 
between  these  two  thickness  measures  as  will  be  shown  subsequently.  Note  that  for  any 
variable  4>,  the  density-weighted  Favre  average  is  denoted  by  <j>  while  the  unweighted 
average  is  denoted  by  (j). 

The  temporal  growth  rate  of  the  vorticity  thickness  in  a  frame  moving  with  the 


convection  velocity  is  defined  by 


_  1  dS» 
Au  dt 


(7) 


and,  for  the  incompressible  case, 


(8) 


where  Cs  is  a  nominal  constant.  As  established  in  the  previous  section,  Cs  decreases 
with  Mc  but  here  the  influence  of  density  ratio  on  Cs  is  of  interest.  In  the  low-speed 
case,  [15],  based  on  limited  low-speed,  variable-density  shear  layer  data,  hypothesized 
that  Cs  is  independent  of  density  ratio.  Our  results  at  Mc  =  0.7  are  shown  in  Table  1. 
It  is  clear  that  the  influence  of  density  ratio  on  Cs  is  insignificant  even  in  the  high-speed 


regime. 

On  the  other  hand,  the  momentum  thickness  growth  rate  decreases  systematically 
with  departures  of  density  ratio  from  its  reference  value  of  s  =  1  as  shown  in  Table  2. 
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Table  1:  The  influence  of  density  ratio  on 
the  growth  rate  of  the  vorticity  thickness 


Table  2:  Density  ratio  effect  on  growth 
rate  ratio. 


s  Cs 

— 

s  DNS,  Sg/5e.i  Model,  Eq.  (9) 

1  0.117 

4  0.137 

8  0.125 

1  1.00  1.00 

4  0.86  0.82 

8  0.59  0.70 

Self-similar  profiles  of  the  Reynolds  shear  stress  are  shown  in  Fig.  5.  The  peak  values 
do  not  vary  with  density  ratio,  consistent  with  the  insignificant  change  in  vorticity 
thickness  growth  rate.  However,  the  location  of  the  peak  Reynolds  shear  stress  shifts 
to  positive  values  of  x2/56  along  with  the  (approximately  coincident)  position  of  the 
dividing  streamline.  The  location  of  the  dividing  streamline,  U  —  0,  can  be  determined 
from  the  velocity  profiles  in  Fig.  6.  Also,  from  the  density  profiles  in  that  figure  it  is 
clear  that  the  location  of  the  dividing  streamline  and  the  peak  Reynolds  shear  stress  lie 
in  the  low-density  side  of  the  shear  layer.  The  shift  of  the  center  of  the  velocity  profile, 
U  -  0,  with  respect  to  the  center  of  the  density  profile,  p  =  (pi  +  p2)/2  =  0.7,  increases 
with  density  ratio  as  shown  by  Fig.  6.  The  DNS  data  are  used  to  parameterize  this 
shift  by  a  quantity  a(s).  This  shift  towards  the  low-density  stream  increases  when  s 
increases  resulting  in  the  systematic  decrease  of  peak  values  of  pR\2  and  therefore  of  the 
momentum  thickness  growth  rate.  An  analysis  has  been  carried  out  to  explain  this  shift 
and  yield  the  following  model, 

i  t~~~t  =  1  -  0.4A(s)a(s)  (9) 

6g(s  =  1) 

where  \{s)  =  (s  -  1)/(1  +  s)  and  a(s)  is  a  parameter  that  determines  the  shift  of  the 
center  of  the  mean  velocity  profile  with  respect  to  that  of  the  mean  density  profile.  Table 
2  shows  that  Eq.  (9)  is  a  good  approximation  to  the  observed  effect  of  density  ratio  on 
the  momentum  thickness  growth  rate. 

2.3  Effect  of  heat  release 

Popular  models  of  flow/chemistry  interaction  such  as  the  unsteady  laminar  flamelet 
[19,  20]  and  the  conditional  moment  method  [21,  22,  23]  are  based  on  the  knowledge  of 
the  statistics  of  two  variables;  a  mixture  fraction  variable,  C,  that  represents  the  state  of 
mixing  and  a  second  variable,  the  scalar  dissipation,  x,  defined  as  2D  VC  •  VC,  with  D 
the  mass  diffusivity.  These  models  attempt  to  simplify  the  coupling  between  the  scales 
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Figure  5:  Self-similar  profiles  of 
Reynolds  shear  stress,  Ru,  for  various 
density  ratios.  The  cross-stream  coordi¬ 
nate.  x2,  is  nondimensionalized  with  the 
momentum  thickness. 


Figure  6:  Self-similar  profiles  of  the 
mean  velocity  and  density.  The  center 
(where  U  =  0)  of  the  mean  velocity  pro¬ 
file  is  shifted  with  respect  to  the  center 
(where  p  =  0.7)  of  the  mean  density  pro¬ 
file. 


of  turbulence  and  the  several  orders  of  magnitude  smaller  chemical  scales.  Since  little 
is  known  on  the  effect  of  heat  release  on  the  mixture  fraction  and  scalar  dissipation,  we 
have  focussed  our  DNS  analysis  on  these  quantities. 

For  obtaining  the  leading-order  effect  of  heat  release  it  is  not  necessary  to  consider 
multi-step,  finite  rate  chemistry.  Thus,  the  single-step,  infinitely  fast  irreversible  com¬ 
bustion  of  methane  with  air  is  considered, 

CH4  +  202  ->  002  +  2 H20  (10) 


The  heat  release  is  represented  by  Ce  =  Taj/To  —  1,  where  Ta(t  is  the  adiabatic  flame 
temperature  and  T0  is  the  reference  freestream  temperature.  A  stoichiometric  mixture 
fraction  value  of  G  =  0.2  was  fixed  by  diluting  the  fuel  stream  with  nitrogen.  In  addition 
to  the  value  of  0e  =  7  corresponding  to  the  chosen  stoichiometry,  a  case  with  0e  =  3.5 
with  the  heat  release  constrained  to  half  the  actual  value  is  also  considered. 

The  probability  density  function  (pdf),  P(G^)j  of  the  scalar,  G  is  a  function  of 
the  normalized  cross-stream  coordinate,  rj  —  in  the  self-similar  shear  layer.  Fig¬ 

ure  7(a)  shows  the  scalar  pdf,  P{ C;t?),  at  7?  =  0.2  (average  position  of  the  flame)  for  the 
turbulent  shear  layer  in  cases  Ce  =  0,  3.5  and  7.  The  values  of  P(C  =  Cs'iV  —  0-2)  in 
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P(C;n=0.2) 


Figure  7:  Scalar  pdf  (a)  Ce  =  0,  3.5  and  7  at  77  =  0.2  and  (b)  77  =  ±0.4.  Here,  77  =  x2/SL 
is  the  normalized  cross-stream  coordinate. 


Figure  8:  Scalar  dissipation  (a)  and  scalar  r.m.s  (b)  for  Ce  =  0,  3.5  and  7. 
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both  cases  with  heat  release  are  about  twice  the  value  for  the  case  without  heat  release. 
Figure  7(b)  shows  the  scalar  pdf  at  r)  =  ±0.4.  The  pdf  is  no  longer  symmetric,  since 
the  stoichiometric  value,  G  =  0.2,  introduces  an  asymmetry  in  the  mixing  layer.  An¬ 
other  feature  of  the  reacting  cases  is  that  external  intermittency  is  suppressed.  These 
changes  should  be  considered  when  developing  models  for  reacting  flows.  Moreover,  the 
use  of  pdf  measurements  that  correspond  to  non-reactive  flows  could  lead  to  incorrect 

predictions. 

Figure  8(a)  shows  the  unconditional  scalar  dissipation,  x,  normalized  by  vorticity 
thickness,  as  a  function  of  the  coordinate  77.  It  can  be  seen  that,  although  the  non¬ 
reacting  and  reacting  cases  have  similar  profiles  when  8U  is  used  as  the  scaling  parameter, 
there  is  some  reduction  in  the  thickness  of  the  unconditional  scalar  dissipation  profile 
when  there  is  heat  release.  Figure  8(b)  shows  profiles  of  r.m.s  (root  mean  square)  scalar 
fluctuations  for  Ce  =  0,  3.5  and  7.  The  levels  of  scalar  r.m.s  with  heat  release  are  lower 
than  the  values  measured  in  the  case  without  heat  release. 

The  pdf  of  scalar  dissipation  conditioned  on  scalar  mass  fraction  has  a  complex 
behavior  in  the  reacting  shear  layer,  a  behavior  which  is  currently  being  analyzed. 

3.  REDUCED  CHEMISTRY  FORMULATIONS 

The  structure  of  hydrogen-air  diffusion  flames  in  conditions  typical  of  high-speed 
propulsion  applications  is  studied  with  detailed  and  reduced  chemistry  mechanisms  [16]. 
Pressures  ranging  from  0.25  atm  to  10  atm,  fuel  stream  temperatures  ranging  from  300 
K  to  COO  K,  and  air  temperatures  from  500  K  to  1500K  are  considered.  It  is  found  that  a 
21 -step  detailed  kinetic  mechanism  can  describe  the  structures,  extinction  and  ignition 
of  hydrogen-air  diffusion  flames  over  this  range  of  conditions.  More  importantly,  the 
efficacy  of  reduced  mechanisms  is  examined  in  this  study.  A  four-step  mechanism  based 
on  steady  state  approximations  for  HO2  and  H2O2  gives  essentially  the  same  results  as 
the  detailed  mechanism.  A  three-step  mechanism,  based  on  an  additional  steady-state 
approximation  for  OH,  and  a  two-step  mechanism,  based  on  additional  steady-state 
approximations  for  both  OH  and  O  is  also  evaluated.  Both  mechanisms,  especially  the 
three-step  variant,  are  found  to  give  acceptable  predictions  of  ignition,  extinction  and 
heat-release,  but  do  not  predict  radical  concentrations. 

A  three-step  reduced  mechanism  for  hydrogen/oxygen  combustion  consists  of  a  fast 

fuel-consumption  step, 

3 i/2  ±  O2  2H2O  +  2 H 
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together  with  two  recombination  steps' 


H  +  H  H2 
2H2  +  02  — ^  2H20  . 

The  above  mechanism  is  valid  for  strain  conditions  ranging  from  weakly  strained  flames 
to  near  extinction.  Using  the  above  mechanism  with  explicit  reaction  rate  terms  in  DNS 
or  LES  is  problematical  because  of  the  fast  fuel-consumption  step.  Therefore,  generalized* 
coupling  functions  have  been  developed  to  give  conservation  equations  that  contain 
explicitly  only  the  slow  radical  recombination  rate  while  the  fast  fuel-consumption  step 
is  treated  with  the  new  coupling  functions  that,  by  virtue  of  being  continuous  functions, 
allow  discrete  differentiation,  without  any  special  techniques,  in  a  simulation.  Species 
diffusivities  that  differ  from  the  thermal  diffusivity  are  permitted.  The  mathematical 
model  corresponds  to  a  reaction  sheet  at  which  fuel  is  consumed  in  a  diffusion-controlled 
fashion,  yet  allowing  finite-radical  recombination.  Good  results  of  the  laminar  flame 
structure  are  obtained  by  [17]  with  this  formulation  when  applied  to  a  hydrogen-air 
laminar  mixing  layer  with  free-stream  temperatures  above  the  crossover  temperature 
corresponding  to  the  second  explosion  limit. 

Generalized  coupling  functions  derived  from  a  three-step  reduced  mechanism  for  the 
combustion  of  methane  [18]  are  currently  being  evaluated  for  application  in  DNS. 

4.  THE  TURBULENT  JET:  DNS 

Turbulent  plane  jets  are  prototypical  free  shear  flows  of  practical  interest  in  propul¬ 
sion  and  other  applications.  While  considerable  experimental  research  has  been  per¬ 
formed  on  planar  jets,  very  few  computational  studies  exist.  Our  work  is  the  first  com¬ 
putational  study  of  a  spatially-developing  turbulent  jet  utilizing  DNS.  Comparisons  of 
two-  and  three-dimensional  simulations  [24,  25,  26]  show  the  inadequacy  of  the  former. 
The  development  of  the  velocity  and  scalar  field  including  successful  comparison  with 
experimental  data  is  discussed  in  [27]  while  the  effect  of  nozzle  conditions  is  reported 
in  [28].  A  Ph.D.  thesis  [29]  has  also  resulted  from  these  studies. 

The  discussion  here  is  limited  to  some  features  of  the  scalar  field  development  in  a 
turbulent  jet  with  initial  Reynolds  number,  Re/,  =  3000,  that  increases  to  Re/,  =  4700 
at  the  outflow.  Here,  Reh  =  Ujh/u  is  based  on  the  uniform  jet  exit  velocity,  Uj,  and  the 
nozzle  height,  h.  Figures  9  and  10  show  instantaneous  isocontours  of  the  passive  scalar 
in  an  xy-plane  (side  section)  and  in  a  rrz-plane  (top  section),  respectively.  The  rrz-plane 
shown  is  at  a  station  y/h  =  0.467  which  is  just  below  the  upper  shear  layer  at  the  inflow. 
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Figure  9:  Instantaneous  passive  scalar  contours  on  an  xy-plane,  z  =  0.0.  The  nozzle 
height  is  denoted  by  h. 


Figure  10: 
shear  laver 


Instantaneous  passive  scalar  contours  on  a  xz- plane  just  below  the  upper 
y/h  =  0.44. 


Figure  11:  Variation  of  the  probabil¬ 
ity  density  function  of  the  passive  scalar 
across  the  jet  at  x/h  =  7.0.  The  cross¬ 
stream  coordinate,  y,  is  normalized  by  the 
scalar  half-width,  5 {.  Symbols:  □  y/^  — 
0.00,  A  y/5i  =  0.43,  <1  y/5^  =  0.87,  0 
yfa  =  1-30,  O  v/h  =  1-74. 


Figure  12:  Variation  of  the  probabil¬ 
ity  density  function  of  the  passive  scalar 
across  the  jet  at  x/h  =  11.5.  Symbols: 
□  y/6z  =  0.00,  A  yfa  =  0.38,  <1  y/5^  = 
0.72,  0  y/5i  =  1.07,  Q  y/5 f  =  1.45. 


At  this  lateral  station,  the  inflow  fluid  is  predominately  jet  fluid.  In  these  figures,  white 
indicates  pure  coflow  fluid,  f  =  0.0,  while  black  indicates  pure  jet  fluid,  £  =  1.0.  The  jet 
column  (almost  pure  jet  fluid)  is  clearly  evident  in  figure  9  and  extends  until  x/h  =  5.0 
where  it  begins  to  break  down.  Shear  layer  and  jet  instabilities  that  have  large  spanwise 
extent  (see  Fig.  10)  are  evident  in  the  initial  region,  0  <  x/h  <  5.  Further  downstream, 
small-scale  structures  with  limited  spanwise  coherence  dominate. 

The  scalar  pdf’s  in  turbulent  jets  have  been  a  subject  of  some  controversy.  Figure  11 
shows  the  scalar  pdf’s  across  the  jet  at  the  station  x/h  —  7.0.  At  this  station,  the  proba¬ 
bility  densitv  functions  for  y/5 ^  —  0.43,  0.87  and  1.30  are  decidedly  non-marching.  While 
the  mean  scalar  values,  solid  shaded  symbols  at  the  top  of  the  plot,  vary  strongly  for 
0.43  <  y/S^  <  1.30,  the  central  peak  in  the  pdf’s  does  not  march,  remaining  stationary 
at  scalar  values  £  ~  0.5.  At  each  of  these  lateral  stations  there  is  a  second  strong  peak 
in  the  probability  density  function  corresponding  to  either  pure  jet  fluid,  y/ =  0.43, 
or  pure  coflow  fluid,  y/5c  =  1.30,  or  both  y/5i  =  0.87.  The  probability  density  functions 
at  y/6(  =  0.0  and  1.74  are  quite  broad  and  show  a  dominant  peak  of  pure  jet  and  pure 
coflow  fluid,  respectively.  The  non-marching  pdf’s  at  this  streamwise  station,  x/h  =  7.0, 
indicates  that,  in  this  region  of  the  jet,  mixing  is  dominated  by  the  large-scale  engulfing 
of  fluid. 

Figure  12  shows  the  probability  density  functions  for  the  passive  scalar  field  further 
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Figure  13:  Evolution  of  the  jet  half-width 
in  LES. 


Figure  14:  Profiles  of  streamwise  turbu¬ 
lence  intensity  in  LES. 


downstream  in  the  jet,  x/h  =  11.5,  after  the  fluctuating  velocity  and  scalar  fields  ap¬ 
proach  self-similarity.  At  this  station,  the  PDF  at  the  jet  centerline  no  longer  shows  a 
strong  peak  corresponding  to  pure  jet  fluid,  £  =  1.0,  implying  that  small-scale  mixing  has 
eliminated  the  presence  of  strong  patches  of  unmixed  jet  fluid.  The  probability  density 
functions  at  y/5^  =  0.0  and  0.38  exhibit  “marching”  type  behavior  at  this  streamwise 
station.  There  is  a  single  peak  in  the  PDF  which  marches  along  with  the  mean  value 
of  the  passive  scalar  at  that  lateral  location.  Near  the  outer  edges  of  the  jet,  there  is  a 
secondary  peak  corresponding  to  pure  coflow  fluid,  suggesting  intermittent,  large-scale 
entrainment. 


5.  THE  TURBULENT  JET;  LES 

Large  eddy  simulations  of  a  spatially  developing  jet  have  been  performed  [30]  and 
three  different  subgrid  stress  models  compared:  the  standard  Smagorinsky  model,  the 
dynamic  Smagorinsky  model  and  the  dynamic  mixed  model.  A  set  of  simulations  at  a 
low  Reynolds  number,  Red  =  3,000  enables  comparisons  with  our  previous  DNS  [27]. 
Then  simulations  are  performed  at  a  higher  Reynolds  number,  Red  =  30,000,  more 
typical  of  engineering  applications.  Fig.  13  shows  the  evolution  of  the  jet  half-width. 
The  standard  Smagorinsky  model  is  much  too  dissipative,  while  the  other  two  models 
compare  well  with  both  DNS  and  experiments.  Fig.  14  shows  that,  in  the  high  Reynolds 
number  jet,  both  the  mixed  and  dynamic  Smagorinsky  models  give  turbulent  intensity 
profiles  consistent  with  experimental  data. 

Our  a  priori  tests  using  the  exact  subgrid  stresses  obtained  by  filtering  the  DNS 
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jet  data,  show  that  the  subgrid  shear  stress  is  severely  underpredicted  by  the  dynamic 
Smagorinsky  model.  On  the  other  hand,  surprisingly,  the  dynamic  Smagorinsky  model 
provides  satisfactory  overall  results  concerning  flow  statistics  in  a  posteriori  calculations, 
as  shown  by  Fig.  13  and  Fig.  14.  Terms  in  the  equation  for  streamwise  velocity,  U i,  are 
examined  to  understand  this  apparent  contradiction.  The  dominant  balance  is  found 
to  be  between  mean  convection  and  the  transverse  gradient  of  the  sum  of  the  resolved 
Reynolds  shear  stress  and  the  unresolved  SGS  stress.  After  examination  of  the  mean 
momentum  balance,  it  is  found  that  the  sum  of  the  resolved  and  subgrid  shear  stress  is 
similar  with  both  models,  thus  explaining  their  similar  performance. 

We  are  currently  evaluating  the  ability  of  LES  with  subgrid  scalar  transport  models  to 
predict  the  evolution  of  mean,  r.m.s.  and  pdf’s  of  the  scalar.  Also,  our  new  developments 
[31]  in  LES  of  nonequilibrium  turbulent  flows  will  be  tested  to  predict  the  well-known 
sensitivity  of  jets  to  different  coflow  speeds,  a  difficult  problem  for  RANS  calculations. 

6.  COUETTE  FLOWS  AND  RELATED  STUDIES 

A  study  of  extinction  by  wall  cooling  and  the  influence  of  compressibility,  heat  trans¬ 
fer  and  shear  on  the  rate  of  creation  of  product  in  laminar  Couette  flow  [32]  has  resulted 
from  the  research  conducted  as  part  of  this  grant.  In  addition  to  being  of  fundamen¬ 
tal  interest,  this  contribution  provides  the  basis  for  development  of  extended  flamelet 
models,  i.e.,  models  applicable  to  high-speed  turbulent  reacting  flows.  One  of  the  fun¬ 
damental  problems  in  the  modeling  of  such  flows  relates  to  the  treatment  of  fluctuating 
pressure  gradients  which,  although  neglected  in  previous  theories,  have  been  found  to  be 
important  in  recent  DNS  studies.  In  [33]  a  revised  theory  of  premixed  turbulent  com¬ 
bustion  involving  a  new  model  for  such  fluctuating  gradients  is  developed  and  applied 
to  turbulent  reactant  streams  impinging  on  a  wall,  a  flow  for  which  a  limited  amount  of 
accurate  data  on  first  and  second  moment  quantities  is  available.  Excellent  agreement 
between  theory  and  experiment  is  achieved.  Particularly  impressive  is  the  agreement 
with  respect  to  the  mean  axial  flux. 

Suggested  by  this  success  is  a  revisiting  of  the  theory  of  turbulent  flames  normal 
to  an  oncoming  reactant  stream.  The  hope  here  is  that,  depending  on  the  intensity 
of  the  turbulence  in  the  stream  relative  to  the  laminar  flame  speed,  a  measure  of  the 
chemical  kinetic  rate  of  the  mixture,  we  shall  predict  either  countergradient  or  gradient 
transport.  It  is  now  known  from  both  DNS  results  and  experiment  that  under  conditions 
of  intense  turbulence  the  countergradient  transport  seen  so  frequently  in  computations 
and  experiment  reverts  to  gradient  form.  This  transition  has  never  been  predicted  by  a 
moment  method,  but  it  is  hoped  that  the  new  model  for  fluctuating  pressure  gradients 
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will  permit  such  a  prediction. 

As  a  follow  on  to  these  studies,  models  for  the  mean  rate  of  creation  of  products 
are  being  currently  examined;  there  are  many  such  models,  but  they  are  never  critically 
assessed,  by  themselves.  However,  in  our  studies  of  flames  in  stagnating  turbulence,  this 
separation  can  be  readily  made,  and  each  chemistry  model  assessed  independently. 
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